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Abstract. We examine the unique spectroscopic features which give rise 
to quantum information processing attributes of one-dimensional J-aggregate 
systems, and as revealed by entanglement measures such as the von 
Neumann entropy, Wootters concurrence and Wei-Goldbart geometric measure of 
entanglement. The effect of dispersion and resonance terms in the exciton-phonon 
interaction are analyzed using Green function formalism and present J-aggregate 
systems as robust channels for large scale energy propagation for a select range 
of parameters. We show that scaling of the third order optical response x^^"^ with 
exciton delocalization size provides an experimentally demonstrable measure of 
quantifying multipartite entanglement in J-aggregates. 
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1. Introduction 

Quantum J-aggregates in organic dye molecular systems reveal an unusually sharp 
and intense red-shifted optical absorption band arising from molecular aggregation. 
This spectroscopic feature provides a spectacular example of nonlocal effects due to 
collective properties of the Frenkel exciton [1-4]. The narrow isolated band of J- 
aggregates usually formed in liquid or solid organic matrix, generally appears adjacent 
to a monomer based broad absorption band incorporating vibrational features. 
The absence of a similar vibrational structure in the polymer J-band is a striking 
phenomenon, almost equivalent to a decoherence-free subsystem, widely discussed in 
quantum information studies [5]. 

It is seen that subtle features of quantum J-aggregates cannot be adequately 
explained using conventional theory of exciton dynamics [6-8]. Recent technical 
advances in spectroscopic techniques such as two-dimensional (2D) coherent electronic 
spectroscopy with fast femtosecond time resolution [9] has made it possible to 
observe intricate dephasing processes of delocalized exciton states in the single J-band 
aggregate system [10]. Such detailed experimental probes of the entangled properties 
highlight the potential in utilizing opto-electronics properties of J-aggregate systems 
for quantum information processing. Optical-based applications have already been 
demonstrated in organic semiconductor microcavities [11]. The potential application 
to innovative "frontier technologies" based on new understanding of the puzzling 
effects of quantum sized structures therefore provide much impetus for further 
examination of J-aggregates, particularly from a quantum information perspective. 

Narrow J-bands are generally studied using a model of identical monomers 
interacting with vibrations [12, 13] or an alternative model based on rigid monomers 
by just varying electronic site energies [14]. The SP parameter, given by the ratio of 
inter-monomer interaction energy to the vibrational width of the monomer spectrum, 
was used by Simpson et al [12] to provide a measure of narrowness of the J-band. In 
the strong coupling region with high SP values, excitation transfer between monomers 
occurs much faster than relaxation due to lattice vibrations. In Knapp's Frenkel 
exciton model [14], specific details of the background medium was discarded with 
attention focusscd on the the J-band. The time of propagation of the Frenkel exciton 
was considered faster than time variations in electronic site energies. The term 
"exchange narrowing" [14] was used to describe the decrease in linewidth by 1/VN in 
the absence of intersite correlation, being the number of monomers. The exchange 
narrowing process can be associated with N monomers participating in a collective 
mode as a multipartite system, each monomer thus experiencing an effective disorder 
of magnitude instead of D, the full disorder experienced by an isolated monomer. 

It is well known that quantum non-locality is a collective property of quantum 
system of two or more subsystem (e.g monomers) in which the dynamics of constituent 
subsystems depend on spatially separated adjacent subsystems [5]. Pauli exclusion 
principle excludes the possibility of two excitons occupying the same molecular 
site, thus processes such as exciton-exciton scattering and excitation/de-excitation 
at a site invariably leads to entanglement in linear excitonic chains. The "spooky 
action at a distance" effect [15] ensures minimal participation of vibrational modes 
in quantum J-aggregates. The time taken for individual monomers to interact and 
form the collective J-band is far less than interaction time with lattice vibration as 
revealed experimentally [9,10]. Recently the puzzling effects of entanglement was 
experimentally demonstrated [16] by joint weak measurement procedures using an 
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entangled photon pair. 

In this work, wc examine the dynamics of J-aggregatcs using a quantum 
information approach in which the origin of its sharp J-band is attributed to 
entanglement and non-local effects associated with a system of excitonic monomers. 
Wc consider that dccohcrcncc effects [17-20] associated with a noisy environment 
invariably disrupts the delicate delocalized excitation which contributes to the 
narrowness of J-aggregate spectrum. In the extreme limit of an infinite chain 
with infinitely strong coupling, wc assume that nonlocal effects is maximal with 
characteristics of an ideal J-band forming an infinitely narrow spike. 

J-aggregates appear remarkably to be the first systems in which entanglement 
effects have been revealed so visibly during the 1930's, understandably Schcibc [2] 
attributed the J-band as arising due to a collective state of electronic excitation. 
Despite the consistent observation of a collective excitonic mode by Scheibe [2], 
Simpson ct al [12] and Knapp [14], attempts to bridge the gap between the well-known 
Frenkel exciton theory and the more recent field of quantum information theory has 
only been achieved recently in our work [21], where we have shown clear and explicit 
links between the quantum search process and actual physical processes taking place in 
the crystal. To simplify the numerical analysis in this work, we utilize the Markovian 
approximation scheme in which the time scales of lattice vibrations are considered to 
be short compared to other dynamical times present in the system. Hence wc neglect 
the flow of information from phonons back to the excitonic system and assume that the 
reservoir system is disengaged from coherence effects between two or more molecular 
sites. In other words, wc consider that the operating times associated with non local 
effects remain unaffected by "memory effects" present in a non-markovian reservoir 
system. We thus distinguish between propagation times of the excitation and times 
associated with entanglement effects of a delocalized exciton which determines the 
duration of coherence between different molecular sites. 

In view of the importance of including entanglement effects in J-aggregates, we 
emphasise the general problem of qiiantifying the level of entanglement in any arbitrary 
multipartite system. Currently there is no satisfactory classification of genuine 
entanglement for a multipartite system of even pure states [22-27]. The manner of 
partitioning individual subsystems affects multipartite entanglement and introduces 
an element of subjectivity when using current classification of entanglement such as 
k-separability. For instance for a selected partition, one can obtain varying degrees 
of multipartite entanglement, changing from a fully separable to fully inseparable 
states with each member of the group sharing genuine multipartite entanglement. The 
possibility that such a change may be non-monotonic was recently shown in an earlier 
work where three members of a group which shared non-zero tangle entanglement 
displayed zero bipartite entanglement between two members [28] . We therefore expect 
investigations related to maximal entanglement of a fragile open systems such as 
excitonic chains with dccohcrcncc effects to remain unresolved for some time. In this 
regard, we point to the possibility that the scaling of the third order optical response 
X^^^ with exciton delocalization size may provide an experimentally demonstrable 
means of categorizing genuine multipartite entanglement, at least in the case of J- 
aggregates. 

This paper is organised as follows. In Sec. 2, we summarize salient features 

of the Frenkel exciton in molecular systems, including a brief description of the 
site basis representation appropriate to excitonic systems. In Sec. 3, we provide 
a brief description of the von Neumann entropy [24,29], and in Sec. 4 we provide 
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mathematical details of derivations based on Green's functions formalism. We also 
present results which show the usefulness of the entropy in capturing the transition 
from coherent to incoherent energy transfer in a finite one-dimensional J-aggregate 
system. In Sec. 5, we use the well-known Wootters concurrence [30] to interpret 
superradiance properties of one-dimensional excitonic systems. Finally, we analyze 
the intricate link between the third order optical response x^^^ a^^d the class of 
entanglement quantity known as the Wei-Goldbart geometric measure of entanglement 
[31] in Sec. 6. 

2. J-aggregates and Frenkel Excitons 

J-aggregates, also called "communal" states [2] , are low dimensional molecular systems 
which display enhanced coherent interactions with light. For instance, a well known 
dye molecule pseudoisocyanine (PIC) forms a quasi-one dimensional arrangement of 
molecules under suitable conditions, displaying the characteristic feature of a sharp 
absorption spectrum red shifted from the broader monomer transition band. It is well 
known that many spectroscopic featiu-es of J-aggregates are interpreted via a model of 
one-dimensional Frenkel exciton system [6], with delocalised excitons considered the 
main reason behind the collective nature of such systems. There have been suggestions 
that the total number of molecules which form such collective states varies from ten 
to several hundreds [14]. 

The optical properties of J-aggregates are invariably linked to geometrical 
attributes of the molecular arrangement, however in this work we restrict our study 
to Frenkel excitons in a quasi one-dimensional system of finite N molecular sites. 
We assume that the one- dimensional Frenkel excitonic model is applicable to one- 
dimensional aggregate systems as well. Excitons in molecular crystals are generally 
treated using the tight-binding model [6] due to the weak overlap of intermolecular 
wavefunctions and hence localization of excitation at individual molecular sites. The 
exact form of these functions are not needed to derive salient properties of Frenkel 
excitons, however every molecular site has an equal probability of being excited by 
an incident photon due translation symmetry, a salient feature which has important 
implications for quantum parallelism. Various mechanisms result in the loss of 
excitonic coherence, in this work we consider the effects due to coupling with lattice 
vibrations as well as the influence of finite temperature. 

The absorption of a photon triggers an initialization process with creation of a 
delocalized exciton, representing an equal weighted superposed form of occupation 
probabilities at all lattice sites. While ideal crystal in vacuum state is equivalent to 
the null state |0)®^ for N sites, photon absorption process initiates a Hadamard- 
like transformation [5] at each site. In the absence of lattice vibrations, the exciton 
eigenvector can thus be written as 

|i^)=Ar-V2^e^i^'i3t|o) (1) 

where the reciprocal lattice vector K = within the Brillouin zone k € [^ y], 
|0) is the vacuum state with all molecules in ground state and Bj is the creation 
operator of exciton at position coordinate I. N is the number of unit cells in the 
crystal which increases with the number of molecules for each unit cell. In Eq.(l), we 
have dropped notations associated with the exciton spin for convenience reasons, and 
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include notations in scalar form as we restrict the model considered here to a finite one- 
dimensional linear chain. It should be noted that we have considered that all molecular 
sites have an equal chance of absorbing a photon to raise it to an excited state |0) 
which can be represented by the Hadamard-like transformation: |0) |0) + |1). It is 
important to note the vital difference between a photon absorption process and a true 
Hadamard gate conversion which can reset a state to its original form. For instance, 
other possible transformations inherent in a true Hadamard gate such as |0) + 1 1) |0) 
and |1) -> |0) — |1) |1) are not present upon successive absorption of two or more 
photons at a molecular site. 

The creator operator of a Frenkel exciton with wavevector K can be easily 
obtained using the Fourier transform of Eq.(l) 

= 7V-V2^e'^-'B/ (2) 

I 

The exciton creation operator b]^ localized in K-sp&ce is delocalized in real space 
[6]. The motion of the exciton in molecular systems is governed by an Hamiltonian 
derived using a tight-binding model [6-8] 

^ex = ^ Ai? + ^ A,m BjBi + Mi^mBjBm (3) 
I \ m^l I m^l 

where AiJ, the on-site (intra-site) excitation energy at equilibrium is the same at all 
sites due to translational symmetry. Di „i is the dispersive interaction matrix element 
which determines the energy difference between a pair of excited electron and hole 
at a molecular site and ground state electrons at neighboring sites [6,7]. Mi^m the 
electron transfer matrix element between molecular sites at I and m. Using Eq.(2), 
we obtain a Hamiltonian diagonal in K space 

iJe. = ^i?o(fc)ijIiJfe (4) 
k 

Eoik) = AE+Y^ Do,m + Mo^rn exp(ifc.m) 

where E{k) is energy of the exciton of wavevector k and subscript denotes the 
absence of lattice site fluctuations. For a simple one-dimensional lattice system, we 
approximate Mi^rn ~ where Vij , the neighbor transfer energy is given by 

Vi,J - ^ (5) 

Hi the molecular transition moment at site i and d is the separation between the two 
sites involved during the energy transfer. 



2.1. Site basis representation 

There are several ways to represent excitonic states as qubit states. For instance, 
Eq.(l) may be considered as the result of a single step of a naturally occurring process 
involving photon absorption, requiring the equivalent of log2 (A'^) + 1 qubits to store 
the superposed state. On the other hand, each molecular site may be occupied by 
an exciton and thus considered a qubit depending on the presence (or absence) of an 
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exciton at a particular site and hence a two-level system. Eq. (3) can be formulated in 
terms of the terms of the exciton density matrix as: 

JV N 

pit) = ^Pn„it)\n){n\ + p„mit)\n){m\ + pl^{t)\m){n\ (6) 

n— 1 n<'m 

where |n) represents a site basis state where the exciton is present at only the 
nth molecular site with all other sites unoccupied. The basis state \n) denotes 
|0)®(n-i) ^ ^ |o)®(A'-n)^ wiicrc N X N is the dimension of the global density 
matrix. |0)(|1)) denotes the absence (presence) of an excitation at a given molecular 
site. The A'' molecular sites occupied by the exciton correspond to the 2^ discrete 
states in the Hilbert space of dimension A^. Transitions between the states in Hilbert 
states hence correlate with the transition amplitudes in the original Hamiltonian in 
Eq.(3). 

A qubit state ^ associated with N > 2 subsystems can be written as generalized 
Schmidt Decomposition belonging to the Hilbert space of He^ in Eq.(3) 

N 

l-Jr) = ^ c„ (|0)®("-i) ® |1) ® |0)®(^-")) (7) 

n=l 

with Schmidt coefficients c^. The Schmidt Decomposition above is well known for 
the separability characterization of pure states. For the simpler ease of just two 
molecular sites N — 2, the bipartite is separable only if there is one non-zero 
Schmidt coefficient. The state can be considered maximally entangled if all the 
Schmidt coefficients are non-zero and equal. The well known Dickc; statc^ \W) = 
(|100...0) -|- |01...0) -I- ... -|- |0...01)) is symmetric under permutations of each basis 
state |n). The absorption of a single photon will give rise to a superposition of all 
possible states as shown in Eq.(l). These collection of states will also include the 
special subset of Dicke states which constitute an equal weighted superposed form of 
basis states. 



3. Quantifying entanglement via the von Neumann entropy 

We use the local density matrix pn ~ Wn|l)nn(^| + (1 — itn)|0)(0| to determine 
entanglement dynamics of the nth qubit associated with a basis which spans the two 
dimensional space (|0) and |1)). The occupation probability at the nth molecular site 
is given by w„ = {K\Bl^Bn\K) where i?„ is the exciton creation operator localized at 
site n and \K) is the exciton eigenvector in if-space (see Eq.(l)). The reduced density 
matrix has a simple form 





A suitable entanglement measure for the Frenkel exciton described by the one- 
dimensional tight-binding model in Eq. (3) and corresponding to the density matrix in 
Eq.(8) is the widely used von Neumann entropy measure [24,29] 

"Sn = -Tr„p„lnp„ (9) 

where 5„ which is zero for pure states, is obtained by tracing over the degrees of 
freedom of the remaining N — 1 molecular sites. The average von Neumann entropy 
measure is then given by 

1 ^ 

5= -5^-Tr„p„lnp„ (10) 

n=l 
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For the ideal delocalized state in Eq.(l), we obtain S = j^lnN at very large N. The 
larger the c^ntropy. greater is the entanglement between excitation at a given site with 
all other remaining sites. On the other hand, 5 = for the localized state and hence 
is a useful measure in characterizing quantum phase transitions at the localization- 
delocalization borders. For the case of an extended state in a finite system of size N, 
we obtain S = -^InN — {1 — j^) ln(l — j^). Using a time dependent density matrix 
(as will be considered in the next Section), we will show that the time evolution of the 
corresponding von Neumann entropy measure can reveal variations in entanglement 
properties as the exciton propagates in a noisy environment. 



4. Green's function formalism and probability propagator 

We begin with the following crystal Hamiltonian [6] 

Ht = Hex + Hp + -ffepl + Hep2 (11) 

Hp =Y.^lUj{q^h, (12) 
q 

k,q 

Hep2 = AT- V2 J2 x(q)BlBk(bl^ + b,) (14) 

where Hex is given in Eq.(4), Hp denotes the phonon energies and b^g{bq) is the creation 
(annihilation) phonon operator with frequency u){q) and wavevector q. Explicit forms 
for the coupling functions F{k,q) and x{q) are given in Ref. [6]. It is important 
to note the difference in the two coupling functions, while x{q) a-cts to localize an 
exciton at its original occupation site, F{k, q) operates during dispersion of an exciton 
to neighboring sites. Thus F{k,q) is influenced by the same factors which determine 
the exciton bandwidth. In situations where xio) >> P{k,q), self-trapped excitonic 
effects are expected to dominate while in the case of comparable coupling strengths 
x{q) ~ F{k,q), coexistence of free and trapped excitons and localization eff'ects are 
expected to prevail. Due to the finite time associated with lattice vibrations, nonlocal 
eflfects which govern entanglement properties of the delocalized exciton are likely to 
experience decoherence effects due to the first coupling term in Eq.(13). At very high 
neighbor transfer energies J, the second dispersive term in Eq.(14) acts as a small 
perturbation on exciton dynamics. We now investigate the effect of these different 
coupling terms using Green's functions formalism introduced by Davydov [6]. 
The Green's function for an exciton at time t is given by [6, 32] 

G{k,t) = ~i{0: ng\T{Bk{t)U{t,0)Bl{0)}\0 : n,) (15) 

Bk{t) =Bk{0)exp{iE{k)t) (16) 

where |0 : rig) denotes a state with Uq as the population of phonons with wavevector q 
and zero exciton population. T is the time ordering operator and U{t, 0) is the Dyson 
function 

U{t,0) =Texp(^-i expiiHrt') Hepi + Hep2 exp{-iHTt')dt'^ (17) 

where Ht, Hepi and Hep2 are given in Eqs. (12), (13) and (14) respectively. A direct 
relation between the Green's function G{k, t) in Eq.(15) and exciton energy E{k) has 
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been derived by Craig et al (see Eqs. 38 and 39 in Ref. [33]) at times larger than the 
inverse of phonon frequencies, 



G{E{k),t) = -iexp 



E{k) + A{k) - N-^Y^ 



t 



X exp[— 7(fc)] exp(— A) 



where 



A(fc)=Ar-i^|F(fc,g)|^ 



k.q 



7(fc) = iV-i ^ \F{k, g)|2 fn, ,5(17,_^) + {n, + 1) ,5(0+^) 



A =N 



huj{q) ± (£(fc + q) 
|X(9)P 



9 



u;{qy 



m) 

{2nq + 1) 



(18) 
(19) 

(20) 

(21) 

(22) 

(23) 



where A(k) yields the principal value of the term associated with it, and the mean 
phonon occupation number at temperature T is given by fiq = {e'xp[huj{q) / ksT] — 
1}~^. We note that A is independent of time unlike the other terms, however this 
only applies in the Markovian approximation. A rigorous treatment in the non- 
Markovian range involves intense numerical computations and therefore will not be 
considered here. Nevertheless the Green's function G(k,t) in Eq.(18) will suffice to 
reveal vital difference in the roles played by the two coupling functions, xio) and 
F{k, q) in influencing the times of coherence between lattice sites. 

The dynamics of a free particle in a one-dimensional system is determined entirely 
by the Green function Gn,m{t) = {n\G{t)\m) [34]. which yields the amplitude for an 
excitation to move from site m at t = to site n at time In a one-dimensional lattice 
with nearest-neighbor interaction energy 2Fcos(ii') and zero lattice vibrations, Gn,o{t) 
is given by 

1 i. r 2vt 

G„,o(t)=Af-'5^e''="exp -i 

k 

Accordingly, the probability P„o(i) that an excitation which originates at site at 
t = and appears at site n at time t can be obtained using Pno{t) = \Gnfl{t)\'^. 
To calculate the probability P„,o(*) fo^' the cxciton in the presence of dispersion and 
resonance terms x{q) and F{k, q) quantifying the exciton-phonon coupling, we use the 
Green's function G{E{k),t) in Eq.(18). For very large N ^ oo, we obtain 



■ cos{K) 



(24) 



2V 



Pnoit) « N-^Jli—t) exp(-A) e 



(25) 



where in{x) is Bessel function of the first kind and ^{k) is given in Eq.(21). In 
the absence of lattice vibrations, Pno(i) reduces to the well-known form obtained by 
Mcrrificld [35] for an excitation propagating on infinite aggregate systems. We have 
assumed that the transfer interaction term V is independent of the environmental 
variables associated with lattice vibrations, also known as Condon approximation in 
Eq.(25). Recent experimental results [9] have shown that clcctron-phonon interactions 
between off-diagonal elements modulate the couplings between the dipoles, resulting 
in decoherence of the transfer interaction term V. Here we assume that the transfer 
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interaction term V remains intact, and Eq.(25) is a reasonably accurate description 

of the cxciton propagation in the presence of lattice vibrations. 

The subtle differences in influence of the coupling functions, x{q) and F{k, q) 
on the probability transfer function P„o(i) is evident from Eq.(25). The dispersive 
coupling factor x{q) tends to localize the excitation at the original site of excitation 
rather than introduce decohcrcnce during propagation. Oh the other hand, the 
resonance coupling term F{k, q) acts to retard the tendency for the exciton to move 
further away from the initial site of propagation by contributing to decohcrcnce effect 
during propagation. Eq.(25) can be rewritten in terms of dimensionless arguments 

P„o(t) « iV-i J2 (ct) exp(-a^) ^ e"*"'* (26) 

k 

where c = , a = 7^7^ and b = ■ Time t is dimensionless as it is scaled by the 
factor u!{q)^^. It is to be noted that t > 1 describes the Markovian limit reasonably 
well. At the Markovian limit, the flow of information from phonons back to the 
excitonic system can be neglected as the exciton hopping times are large compared 
to the phonon reservoir correlation times. This means that two propagating exciton 
states cannot be distinguished as they remain unaffected by correlations within the 
phonon reservoir at times larger than phonon correlation times. 

In order to obtain an explicit result from based on Eq.(25), we make a few 
approximations, firstly we neglect the dispersion in phonon frequencies so that a 
single frequency, uj{q) ^ 50cm~^, representative of the lowest vibrational branch 
in common molecular crystals [33, 36] is considered. The band averages of the two 
coupling functions, xil) and F(k,q) are then varied via parameters a and b (given 
below Eq.(26)), together with the electronic transfer energy, quantified by c. As 
mentioned earlier, we consider a system of a finite number of molecular sites N , and 
evaluate the von Neumann entropy measure of the system using 

N 

5= ^-Tr„p„ln/9„ (27) 

n=l 

The influence of the coupling functions obtained using Eq. (26) on the von Neumann 
entropy measure given in Eq.(27) can thus be easily evaluated. 

In Fig. la, we show results of S versus time t obtained using Eq.(27) for various 
values of N, and for an initial state, Pno(O) = 5n,o- At time t =0, average entropy 
5 =0 due to presence of excitation at just one site, however the entropy rises rapidly 
reflecting the spread of excitation to more sites with constant velocity (determined 
by c). Infact S appears to display an approximate linear dependence on A'^, as it 
should be as a system with large N results in the occupation of more lattice sites and 
greater chance of existing as an extended system. However as we arc considering a 
finite system, the propagating exciton leaves the system under study after a critical 
time. Hence the entropy S can be seen to drop in all three cases of N considered 
in Fig. la. This drop is most pronounced at small N where there faster exit of the 
exciton. It is to be noted that both parameters a and b are set to zero in Fig. la so 
that phonon-coupling processes are absent, and this confirms that the noticeable drop 
in S is indeed due to the finite size effects of our chosen model. 

It is important to note that at small time t, there is a limit to the number of 
sites visited by the exciton and thus we expect an upper limit to S which is noticeable 
at t 2. To illustrate this point, we show results of S versus number of sites N at 
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times t =2 (solid line), t =5 (dashed) and t =9 (dotted) in Fig. lb. We note that at 
small t, entropy S becomes flat after propagation to a finite number of sites. Thus the 
entropy S increases with A'' as more sites are visited before becoming independent of 
N > N* . The critical N* is is dependent on c and hence the transfer interaction term 
V. At larger t, more sites have a non-zero probability of occupation and therefore N* 
increases with t. The dependence of N* on speed of excitation transfer (quantified by 
c) is illustrated in Fig. Ic where we note increase in S and N* with increasing c. The 
effect of the phonon coupling function b = ^^^^^ on S is illustrated in Fig. Id, which 
shows the gradual decrease of S due to increase coupling with phonons. Figs, la-d 
thus illustrate clearly the interplay of three competing effects: 1) the increase of S 
due to spreading; 2) subsequent decrease of S due to finite size effects (propagation 
out of the system) and, 3) decrease of S due to coupling with phonons. 

In Fig. Id, we have included a dotted-dashed line which corresponds to a fully 
extended state at which S = jj In A^— (1 — -i-) ln(l— -i). The entropy S of the extended 
state is higher than other S values calculated for the case of the one-dimensional system 
considered in Figs, la to d. The results indicate that at large N, the entanglement 
persist for longer times and 5* is almost comparable to that of a fully extcndcid system. 
This suggests the large scale energy transport abilities of one-dimensional excitonic 
systems. We next examine in detail the mechanisms associated with varying coupling 
parameters a and b next ( Fig. 2a,b,c.) 

Results obtained by varying a, b, c is illustrated in Fig. 2a,b,c and imply agreement 
with localization tendencies associated with coupling parameters, a and b. The fall in 
S is noticeably steep for high values of 6 ~ 1 at intermediate N ^ 100. This shows the 
higher localization properties associated with the coupling function F{k, q) compared 
to x(g) (see Fig. 2b). Fig. 2c illustrates that high electronic transfer energies or c lead 
to fast rise in S in the time domain (0 < i < 6), which however is not sustained as 
excitation leaves the system as quickly as it spreads. For lower c, the excitation remains 
extended for comparatively longer times with noticeable oscillations associated with 
"death and birth" type excitation and de-excitation process [37] at a monomer site. 
The results obtained in Figs. 1 and 2a,b,c highlight the prominence of the J-band in 
large systems and show the usefulness of the entropy in capturing the transition from 
coherent to incoherent energy transfer. 

So far we have only considered the influence of two coupling terms on the von 
Neumann entropy measure using a simplified model, in which a single excitation is 
considered to propagate to neighboring sites. In reality, the electronic excitation is 
not confined to a single site but distributed over all A'" molecular sites of the crystal. 
The interaction between these sites then gives rise to N states, with each TV state 
characterized by a wavevector k and associated energy level. The superposition of 
several such states give rise to an excitation which propagates from one site another. 
This highlights the inherent difficulties associated with a full computation of von 
Neumann entropy of a realistic system. One has to incorporate the possibility that 
each site becomes a source of excitation that propagates, while receiving an amplitude 
(of the form in Eq.(24)) from other sites. A full superposition of amplitudes from all 
possible sites have to be then taken into account with time as an added dimension to 
obtain an accurate picture of the excitation dynamics. In order to simplify analysis 
of the full model where all sites have equal probability of being excited, we consider 
the Wootters concurrence [30] which provides a convenient means to study bipartite 
entanglement in the next Section. 
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Figure 1. a) The von Neumann entropy S versus time t at N = 200 (solid line), 
N = 100 (dashed) and = 50 (dotted) at a = 6 = and c = 30. 

b) The von Neumann entropy S versus time AT at t = 2 (solid line) , t = 5 (dashed) 
and t = 9 (dotted) at a = 6 = and c = 30. 

c) The von Neumann entropy S versus time at c = 10 (solid line), c = 20 

(dashed) and c = 40 (dotted) at a = = and t =2. 

d) The von Neumann entropy S versus time A' at 6 = (solid line), b = 0.3 
(dashed) and b = 0.5 (dotted) at a = 0, c = 30 and t = 6. Dotted-dashed line 
corresponds to a fully extended state for which S = -^\nN — (1 — ln(l — -^) 



5. Molecular cooperativity as a paradigm of Collective entanglement 

A notable optical response exhibited by one-dimensional excitonic systems is their 
enhanced rate of spontaneous emission compared to that of the monomer, an 
effect known as supcrradiance [38 -40]. In a series of works [41-44], Mukamel and 
coworkers have shown a connection between cooperative emission mechanisms such 
as superradiance and "molecular cooperativity" , a term used to describe the effective 
coherent dynamics of a system of dclocalizcd Frcnkcl cxcitons. The concept of a 
group of molecules acting in tandem to produce enhanced features in x^^\ the third 
order optical response of one dimensional molecular systems was also shown in an 
earlier work by Hanamura and coworkers [39]. The oscillator strength of the K = 
exciton state is directly proportional to the number of molecular sites in the absence of 
decoherence mechanisms, which implicates a scenario where all molecular dipoles act 
collectively in phase to produce giant oscillator strengths. These results highlight the 
potential in harnessing x^"^) as a quantitative measure of multipartite entanglement 
if a quantum-information perspective is attached to the observed "cooperative" 
phenomena. This connection will be further examined in Sec. 6. In this section 
however, we utilize a popular entanglement measure, Wootters concurrence [30] to 
interpret superradiance of one-dimensional excitonic systems. 
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Figure 2. a) The von Neumann entropy S versus time t at N = 150, b = and 
c =10. The four curves correspond to (from top to bottom) a = 0, 0.3, 0.7 and 
1.5. 

b) The von Neumann entropy S versus time i at TV = 100, a = and c =20. The 
three curves correspond to (from top to bottom) 6 = 0, 0.5 and 1. 

c) The von Neumann entropy S versus time t at N = 200, a = 0.5 and b = 0.3. 
The three curves correspond to c = 40 (dotted), 20 (dashed) and 5 (solid line). 
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The well-known Wootters concurrence [30] C provides a measure of entanglement 

between a pah of qubits, it is zero for separable states and is equal to one for maximally 
entangled states such as the Bell states. C has a simple form The average pairwise 
concurrence (C) [45,46] associated with a system of single-particle states is given by 

C is a measure of delocalization and is also known as inverse participation ratio (IPR) 
[47]. Small values of ( correspond to localized states while larger values are linked with 
delocalization. At the extreme limit of localization at a single site C = 1; (C') = 0, 
while for states which are delocalized completely, ( = N which yields the maximum 
entanglement of j^. For extended excitonic system of J-aggregate structures, we 
consider C, as playing similar role as Nc , the characteristic coherence size. A^c has 
been formulated in terms of the density matrix as [44] 

-1 / \ 2 

|2 



Nr.= 



mn\ 




(29) 



where N is the total number of molecules or monomers and pmn is the reduced 
density matrix associated with sites m and n. Nc provides an approximate measure 
by which exciton at site n is entangled with site m, and decays monotonically as a 
function of (n — m). In the absence of any entanglement with extreme localization, 
Pmn = N^^Snm, Nc ~ 1 and thc scalcd concurrence obtained by dividing Eq.(28) with 
jr , (C')s = 0. In the limit of complete delocalization, pmn — 

iV-\ Nc = N which 

yields the maximum (C)s = 1. Thus the inverse participation ratio ( corresponds well 
with the characteristic cohcircmcc; size; A'',.. 

Using a coarse grained approximation (CGA) (described in Appendix B of 
Ref. [43]) to solve the exciton Green's function (similar to Eq.(18)), Spano et al 
[43] determined thc coherence size Nc using a microscopic definition G{0,t) ~ 
exp(— A^ 7rt). The radiative rate is given by A^ in the superradiant limit at which 
excitation associated with individual sites act collectively on the molecular aggregate. 
In particular, Spano ct al [43] obtained an empirical relation relating Nc and thc 
coupling function F{k,q), temperature T, nearest neighbor transfer energy V and 
phonon energy, hjj{q). By rewriting the empirical relation obtained by Spano et al 
[43] in terms of the dimensionless constants b and c (given below Eq.(26)), we obtain 

where the dimensionless temperature term, tu = ^^|y- The molecular cooperativity 
due to the monomer's dipoles lining up collectively in phase, results in a large net 
amplitude, giving rise to superradiant decay [41] and giant nonlinear optical properties 
[40] . The coherence size of the associated J-band is generally taken as the ratio of its 
radiative decay rate to that of a single monomer, however off-diagonal elements of the 
exciton density matrix in Eq.(6) can be also used as basis to compute the coherent size. 
Results obtained by Spano et al [43] showed pronounced decrease in A^c at resonant 
points given by 

A^c ~ 47rc resonance points (31) 

when acoustic phonon frequencies match those between thc ATth exciton level and 
K = Q exciton level. The significance of these resonant points is attached to a transfer 
of oscillator strength away from the K = Q state due to exciton-phonon scattering. 
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The results obtained using Eq.(30) in terms of the concurrence measure given in 
Eq.(28) is illustrated in Fig. 3. The figure shows gradual decrease in concurrence C 
with aggregate size N for various values of b, c at constant t/. = 2. The temperature 
dependence of Nc in Eq. (30) accounts for the sensitivity of lifetimes of the J-aggregate 
at increasing temperatures. Increasing the coupling term F{k, q) (or b) leads to 
destruction of intermolecular cooperativity, hence there is associated decrease in 
concurrence C as shown in Fig. 3. The superradiance phenomenon is thus less likely to 
occur at large b values. Wc note that higher concurrence C is obtained in systems with 
faster propagation of excitation (higher c values) as also shown in Fig. 3. It is imporant 
to note that we have restricted consideration to only the pair-wise concurrence C here, 
and not multipartite entanglement, which will be considered in the next section. 
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0.6 
0.4 
0.2 



0.0 

60 80 100 120 140 160 180 200 

N 

Figure 3. Concurrence C versus aggregate size N at ti^ = 2. Solid curves for 
c = 15, correspond to (from top to bottom) b = 0.5,0.1. Dashed curves for c = 
5, correspond to (from top to bottom) b = 0.5,0.1. 



6. Theory of Two-exciton states 

In this section, we examine the optical features of J-aggrcgatcs which present an 
experimentally quantifiable measure of multipartite entanglement, in the presence of 
dcphasing mechanisms [20]. In this context, the class of entanglement quantity known 
as the Wei-Goldbart geometric measure of entanglement [31] may be convenient in 
analysing multipartite entanglement of J-aggregates based on properties of the third 
order optical susceptibility, x^"^-*. An important feature related to non-linear optical 
properties of material systems is the existence of one- and two-exciton states. One- 
exciton states are easily identified as states in which an excitation is present at any 
one of the N monomers while in a two-exciton state, any two of iV monomers are 
excited while the other monomers remain in the ground state. We emphasize that the 
two-exciton state here is vastly different from the biexciton state which consists of two 
coupled exciton pair system. 

Two-exciton states may be identified with entangled states formed due to 
superposition of two spatially uncoupled excitons, such states having twice the 
bandwidth of the usually considcircid onc-cxciton states. This can be shown by writing 
the diagonal Hamiltonian of the two-exciton state in K space (without involvement 
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of phonon operators) as 

ffex = ^Ei{k)Bl,^Bi^k + E2{k)Bl^B2,k + Tfc(zj) (32) 

k 

where Ei{k) is energy of the exciton (labelled by subscript 1) with wavevector k and 
E2{k) is energy of the second exciton (labelled by subscript 2) with wavevector k. 
The exciton operators B^ j,, Bi j~ retain the same meaning as used in earlier sections. 
Vmn denotes the interaction energy for coupling of transition moments for molecules 
at m,n. Eq.(32) can be diagonalized further in terms of the exciton operators 

Xa{k) = cosp{k)Bl^ + smp{k)Bl^ (33) 
Xb{k) = - sin/3(fc)Bi,fe + cos/3(A;)B2,fc 

we obtain 

^e. = J^E:{k)xa{k)^Xa{k) + El{k)xb{k)^ Xb{k) (34) 

k 

where i?*(0) and -E^(O) are exciton energies associated with Davydov branches at 
fc = 0. The asterisk in the superscript denotes energies which arc split due to the 
intermolecular interaction. Eq.(34) shows that two-cxciton states arc highly entangled 
with the degree of entanglement determined by cosP{k) and sin /?(fc). We next make a 
connection between two-exciton states and generalized symmetric states via the third 
order optical susceptibility, x^"^^ • 

6.1. Third order optical susceptibility, x^^^ o-nd entropy measure 
The third order optical susceptibility, x^^^ can be written as [40] 

X^'^\-l^s;i^\ -,1^1,1^3) = X! \x'-c\-'^s\'^l,'^2,'^i)+x'i*\-'^s\'^l,'^2,(^i) (35) 

where the summation is performed over all possible permutations of frequencies of 
incident light, Wi, a;2, and = Ylii'^i- Expressions for Xc(—Ws;cji, 072,^3) 
and its conjugate are lengthy and given in an earlier work by Knoester [40], in 
terms of ^o;ei and /Xe2;e2,e3- Mo;ei IS the dipole moment associated with transition 
between the ground state and one-exciton state in which an excitation is present at 
any one of the N monomers, a generalization of the well known Dicke state \W) = 
(|100...0) + |01...0) + ... + |0...01)). ^e2;e2,e3 is the dipolc moment associated with 
transition between the one-exciton state and the two-exciton state. Using Eq.(34) as 
basis, we consider that the two-exciton state are a manifestation of symmetric states, 
which are dependent on the number of ground states or O's [31,48] 

\S[N,M)) = sl^^^^^ |0)®^®|1)®(^-^) (36) 

permutations 

The states in Eq.(34) are known to have an entanglement entropy measure [31,48] 

Ea{N,M)= -\nK{N,Mf (37) 



=V m!(7V-M)! U) • 



Quantum information processing attributes of J- aggregates 



16 



For fixed N, the minimum entanglement eigenvalue or geometric measure of 
entanglement [31] A^ax occurs for M = N/2 (even N) and M = {N ± l)/2 (odd 
N). 

We consider that transitions involved in nonlinear optical effects occur via a 
simplified pathway represented as follows 

\S{N, 0)) — ^ \S{N, 1)) (one - exciton state) (39) 
— > |S'(7V, 2)) (two -exciton state) 

The terms of the form given in Eq.(36) are expected to give rise to the observed 
dependence of x^^' on non-local effects associated with energy transfer within various 
exciton site basis. Accordingly x'^^^ scales with the; entanglement entropy E\ in 
Eq.(37), however the explicit nature of relation between the two quantities is not 
immediately clear, this predicament also applies to a possible relation between x^^' 
and genuine multipartite entanglement involving all N monomers. Here we keep the 
analysis tractable by examining the dependence of the intensity of third-harmonic 
generation (THG) Ithg oc \x''^\—'ii^;i^,i^,i^)\'^ on A'' instead. We simplify the 
approach via the following proposed relations 

MO;el ~ ^/""Cl (40) 
/"e2;e2,e3 ~ ^^^^(2 

^ ^A(Ar.l) 
Ea{N,N/2) 

Ej,{N,2) 
Ej,{N,N/2) 

where /i is the transition dipole of individual monomers and the entanglement entropy 
measure E\{n,k) is given in Eq.(37). (i and (,2 are entropy measures associated with 
the transition efficiencies of one-exciton state and two-exciton state respectively. 

Fig. 4 shows gradual decrease in these entropy measures with increasing aggregate 
size A'', with (,2 the entropy measure associated with two-exciton band seen as more 
resilient, which is expected as the associated entanglement persists for longer times 
due its wider bandwidth. In both cases, the entropy decreases with N due to the faster 
increase of Eji^{N,N/2) with A^. Using Eq.(40) and w = |Ai? where the excitation 
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Figure 4. Entropy measures C^i (solid line) and 1^2 (dotted) versus aggregate size 
N. 
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energy AE is defined in Eq.(3), we write the absolute value of the third order optical 
susceptibility with damping rate F as 



(41) 



2Th^ (w2 - 

Fig. 5 shows gradual increase of ;^|x^"^H~3w; w, w, w)| (in which a constant term in 
Eq.(41) was factored out) with aggregate size N before reaching an almost constant 
value which is independent of the size of system. The results of Fig. 5 illustrate 
the important fact that properties of the J-band is determined by a critical number 
of coherently coupled and entangled monomers rather than the total number N of 
monomers present in the system. 

The slope in Fig. 5 yields the degree of entanglement or multipartite entanglement 
measure. In this respect, the onset of zero slope during change of third order optical 
susceptibility with the system size TV can be viewed as the point of maximum entangled 
state. The results in Fig. 5 are consistent with earlier results [40] which showed the 
critical role played by the exciton delocalization length even Knoester [40] did not 
analyse his results from a quantum entanglement perspective. 

The results obtained in Fig. 5 show that scaling of the third order optical 
response x'^^ with exciton delocalization size may be used as basis in quantifying 
multipartite entanglement in J-aggregatcs or systems composed of excitonic spin 
chains. The result obtained here is significant considering that a rigorous definition of 
genuine multipartite entanglement is still lacking in the literature. Thus entanglement 
measures may be defined based on experimental observations and taking into account 
decoherence mechanisms [17,18], exciton-photon interactions and laser detuning 
processes. In this regard, the results obtained so far has established some basic 
guidelines for an experimentally demonstrable measure of identifying multipartite 
entanglement in excitonic systems. In the table below, we list various quantum 
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Figure 5. -^\x^^^ {—3u>;u>,u>,oj)\ ( see Eq.(41)) versus aggregate size N. 
information processing attributes in J-aggregates systems. 



Table 1 Quantum information processing attributes in J-aggregates 
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State/Event 


Quantum information attributes 


Onc-cxciton state 


Encompasses Dickc State 


Two-exciton states 


Encompasses Symmetric States 


Molecular cooperativity 


Paradigm of collective entanglement 


J-aggregate systems 


Channels of quantum information/energy propagation 


Third order optical susceptibility 


Viable multipartite entanglement measure 



In conclusion, we have demonstrated the usefulness of the von Neumann entropy, 
Woottcrs concurrence and geometric measure of entanglement in characterizing the 
quantum information processing attributes quantum J-aggregates. We identify two- 
exciton states (entangled states of uncoupled excitons) and molecular cooperativity 
(manifestation of multipartite entanglement) as features which can be exploited for 
information processing. The entanglement and quantum correlations associated with 
a delocalized exciton over many monomer sites provides a clearer understanding 
of the processes behind observed behavior of J-aggregates. The delocalization of 
electronic states which give rise to narrowness of J-aggregate spectrum as expected 
is complicated by the fine interplay between excitation transfer mechanism and 
decoherence associated with a background of phonon bath. With control of 
the dispersive coupling {x{q)) and resonance coupling (F{k,q)) parameters via 
temperature and matrix characteristics, J-aggregate systems may be utilized as robust 
channels for large scale energy propagation. 

We have also demonstrated that scaling of the third order optical response x^^^ 
with exciton delocalization size provides an experimentally demonstrable measure of 
quantifying multipartite entanglement in Frenkel excitonic systems, x^"^) may be used 
to distinguish between system with low and high fidelities. In future works, we seek 
to address the question of whether the J-band is a maximally entangled state or one 
that exhibits genuine multipartite. Lastly, the results of this work will be useful 
in understanding many-body correlated dynamics of excitonic systems, with impact 
on future experimental work involving measurement of multipartite entanglement. 
This is of considerable significance as a rigorous definition of genuine multipartite 
entanglement still remains an open question. 
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